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Abstract: We derive expressions for strangeness-changing scalar form factors, which 
incorporate known theoretical constraints both at low and high energies. Their leading 
behaviour in the limit of a large number of colours is calculated from the resonance chiral 
Lagrangian. A set of short- distance constraints on the scalar resonance couplings is ob- 
tained, imposing the form factors to vanish at infinite momentum transfer. Making use of 
previous results on S-wave Ktt scattering U, and a coupled-channel dispersion-relation 
analysis, the K7r, Kr] and Krj' scalar form factors are calculated up to 2 GeV. These scalar 
form factors are a key ingredient in the determination of the corresponding scalar spectral 
function which is important in the extraction of the mass of the strange quark from QCD 
sum rules as well as hadronic r decays. 
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1 Introduction 



QCD currents are a basic ingredient of the electroweak interactions. A good understanding 
of their associated hadronic matrix elements is required to control the important interplay 
of QCD in electromagnetic and weak transitions. In this work we shall investigate a simple, 
although highly non-trivial example of such matrix elements, namely strangeness-changing 
scalar form factors. 

Chiral Perturbation Theory (xPT) 0-0 provides a very powerful framework to study 
the low-energy dynamics of the lightest pseudoscalar octet. The chiral symmetry con- 
straints are stringent enough to determine the hadronic matrix elements of the light quark 
currents at very low energies. Unfortunately, these chiral low-energy theorems only apply 
near threshold. To describe the resonance region around 1 GeV, additional dynamical 
information is required. 

One can construct a chiral-symmetric Effective Field Theory with resonance fields 
as explicit degrees of freedom H3- Although not as predictive as the standard chiral 
Lagrangian for the pseudo-Goldstone mesons, the resonance chiral Lagrangian || turns out 
to provide interesting results, once additional short-distance dynamical QCD constraints 
are taken into account ||. At tree level, this Lagrangian encodes the large- iV c properties 



of QCD |[L0|, pT|, in the approximation |L2| of keeping only the dominant lowest-mass 
resonance multiplets. 

In the following, we present a detailed analysis of the scalar strangeness-changing form 
factors, within the resonance chiral framework. These form factors are needed to improve 
the strange quark mass determination from QCD sum rules for the corresponding scalar 
current [|l3| -|20lj . The AS = 1 scalar form factors also govern the leading J = contributions 
to the Cabibbo-suppressed hadronic r decay width. A good theoretical understanding of 
these quantities could result in a more precise determination of m s from r decays |21] ^6 . 

In order to work with a renormalisation-group invariant object, we define the scalar 
form factors f x (s) from the divergence of the corresponding vector current matrix ele- 
ments, 

(0\d»(s lfl u){0)\X) = i(m s -m u ){0\(su)(0)\X) = -i^C x f x {s), (1.1) 

where Axn = M\ — M% and s = p x . Moreover, we have factored out the normalisation 
constants C x , so that f x (s) = 1 at lowest order in x?T. Then the C x for the lowest lying 
hadronic systems K + tt°, K°ir + , K + i] s and K + r]i which contribute to f x {s) are given by 

1 [2 
Ck+tt° = 1, Ck°tt+ = v2, C K + m = — -j=, C K + m = 2y-. (1.2) 
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Below, we shall only work in the isospin limit with fxn = fx+n — fi< n+- 

The form factors associated with the physical 77 and rf fields are easily obtained, taking 
into account the mixing between the two isoscalar states: 

77 \ / cos 9 — sin 9 
rj' J \ sin 6* cos 9 

We shall adopt the normalisation C K + V = C K + m and C K + V i = C K + rjl . Therefore, 

fK V (s) = cos 9 f Km (s) + 2V2 sin 8 f Krn (s) , 

f Kv '(s) = cos9f Km (s) - -^j= sin9f Km (s) . (1.4) 

For our numerical analyses we shall use sin# = —1/3 ~ —20°, which is in the ball-park 
of present phenomenological determinations. Notice that for this particular value of the 
mixing angle /^(s) = to lowest order in x?T. Although this cancellation is no longer 
true at higher orders in the chiral expansion, it indicates a strong suppression of the Krj 
current matrix element. 

We will first derive the theoretical predictions for fx(s) in the limit of a large number 
of colours. This can be achieved in the framework of the resonance chiral Lagrangian. 
Imposing these form factors to vanish at infinite momentum transfer, we obtain a set of 
short-distance constraints on the chiral couplings of the scalar meson resonances. 

The ;\TT loops are subleading corrections in the 1/N C counting. They incorporate the 
unitarity field theory constraints, in a perturbative way, order by order in the momentum 
expansion. For higher energies around the resonance region and above, the low-energy 
expansion breaks down because the unitarity corrections due to re-scattering are no longer 
perturbative. Therefore some kind of resummation of those chiral loops is required to 
satisfy unitarity [p7|-|38|. This effect appears to be crucial for a correct understanding of 
the scalar sector, because the S-wave re-scattering of two pseudoscalars is very strong. 

In a previous paper [[I]], we have presented a detailed study of S-wave Kn scattering up 
to 2 GeV, within the same resonance chiral framework supplemented with the unitarisation 



procedure developed in refs. p3|,p9|1 . We will use those results to perform a calculation of 
the scalar form factors from dispersion relations, explicitly obeying unitarity, also taking 
into account coupled-channel effects. 

In section 2 first the scalar form factors are calculated at tree level in the %PT framework 
including resonances. Then their expressions in conventional at the next-to-leading 
order are reviewed and compared to the resonance %PT approach. In section 3 we derive 
constraints on the scalar form factors by exploiting unitarity and analyticity. These con- 
straints result in a set of coupled dispersion integral equations for the form factors. The 
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Omnes solution of the dispersion relation in the elastic single channel case is also discussed. 
Section 4 contains our numerical analysis for the single as well as coupled-channel cases 
and our central results for the scalar form factors are presented. Finally we close with 
some conclusions and an outlook in section 5. 



2 Effective Lagrangian Results 



2.1 Resonance Chiral Lagrangian 

In the limit of an infinite number of quark colours, QCD reduces to a theory of tree- 
level resonance exchanges fllO|, . At low energies, the dominant effects are governed 
by the lowest-mass meson multiplets, and can be analysed within the resonance chiral 
Lagrangian framework developed in refs. [§, |9|. We refer to those references for details 
on the chiral formalism and notations. We will use a 11(3)^x11(3)^ effective theory with 
nine pseudoscalar Goldstone bosons, which is the appropriate framework in the limit N c 
to infinity ||-||. 

The tree-level calculation of the relevant scalar form factors is straightforward. One 
obtains: 
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1 + 
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^C m (c m Cd) 
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where q and c m are the couplings of the scalar multiplet to the Goldstone bosons, in the 
lowest-order chiral resonance Lagrangian of refs. f|0, Ms is the scalar resonance mass and 
/ m f n = 92. A MeV is the pion decay constant. The iso-singlet momentum dependence p 2 m 
refers to the appropriate meson mass squared in the linear combinations ( |1.4| ), associated 
with the physical rj and rj' fields. 

The numerical values of the couplings q and c m are not very well known. We can get 
a theoretical constraint by enforcing the scalar form factors to vanish at large momentum 
transfer. This seems a very reasonable phenomenological assumption for a composite 
object. All three form factors are zero when s goes to infinity, if the following two conditions 
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are satisfied: 

4c rf c m = f 2 , c m - c d = . (2.2) 

This implies 

c d = c m = f/2 « 46MeV, (2.3) 
and a dipole structure for the scalar form factors, 

= (2 ' 4) 
Taking sin# = -1/3 and M 5 = M K * « 1.4 GeV, 

/^(O) = 1, f Kv (0) = 2V2^ « 0.33, 

3(M& -M„ 2 ,) + A^ 



M'(0) 



2^ 



3 



^0 



0.74. (2.5) 



Notice that /^(0) goes to zero in the limit Mr- = M^. 

The above results can easily be generalised to take into account the exchange of iV 
different scalar multiplets with parameters Ms,i, Cd,i and c m ^ {% = 1, • • • , N). In that case, 
fx( s ) — 1 contains a sum of iV contributions like the ones in eqs. Q2.1|), with the appropriate 
changes on the scalar parameters. The short-distance requirement that fx{s) goes to zero 
for s — > oo then implies the constraints: 



N N 



4^ c <M c «v = f 2 , TJT (°m,i ~ c<i,i) = 0. (2.6) 



, 1 M s 

1=1 1=1 o ; 



The consequences of these conditions have already been investigated in the analysis of S- 
wave Kir scattering, performed in ref. [Q. A more detailed discussion of the phenomeno- 
logical implications of the relations (2~B) in the case of two scalar resonances will be given 



at the end of the next section. 



2.2 Chiral Perturbation Theory 



The Kir and Krj 8 scalar form factors have been computed at the one-loop level in SU(3)lX 
SU(3)# xPT [Q. The result can be written as: 



f 2 8/ 
1 _ A 2 \ — . . s 



f Kn {s) = l + —^± s + —[5s-2E Kw -3-^)J K7r {s 



4I5(ai)_ 1 f K _ ^ Q AL 
8, 

3s - 2T, K7r - 1 Jk V s(s) + -r-f — - 2/i K - 3/j, m ) , (2.7) 



24/2 i — « s i — »w 4A 
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<~ M k j\ , 4L5(a*)„ , 3 / ^ Ai 




f Km {s) = 3 ~~r <; 1 + + -- I 3s - 2S^ - ) J^( 3 



where = M|- + M%. 

The unitarity corrections associated with the re-scattering of the final pseudoscalar 
particles are incorporated through the loop functions Jpq(s). Their explicit expressions 
are given in appendix A. Although the contributions from Goldstone loops are next-to- 
leading in the 1/N C expansion, they generate large logarithmic corrections to the scalar 
form factors. The loops introduce a dependence on the renormalisation scale /x in the 
renormalised chiral coupling L r 5 (fi) and through the explicit factors 



Hp 



M 2 P 



32vr 2 / 2 



In (Mp/jj 2 ) . (2.9) 



Nevertheless, the scalar form factors are of course independent of /x. 

In the large- N c limit, the expressions for the form factors /^(s) and /k>, 8 (s) reduce 

to 

= 3(^lW "* (,) = 1 + 7^ (2 ' 10) 

The comparison^ with the resonance-exchange results in eqs. (|2.1| ) gives the corresponding 
scalar contribution to the 0(p ) chiral coupling L5: 

T s _ . c mCd _J_ (0 ^^\ 



The resonance propagators appearing in (2T) provide an explicit resummation of local 



terms to all orders in the chiral expansion, improving the result ( |2.10| ) to the dipole form 
(2.4). Adding a second scalar resonance with mass Mj, and couplings c' d and c' m , and also 



taking into consideration the 0(p 4 ) coupling constant L 8 , the contributions from the scalar 
resonances take the form 0: 

S = CjPrn ^ ^5 = 4 + C 

5 Mj Mj, 8 2M| 2M|, V ; 

Now, we are in a position to discuss phenomenological consequences of the above re- 
lations for the scalar couplings. Using the second of the short-distance constraints (|2.6| ), 



1 The Kijs scalar form factor in ( |2.l|) includes an additional term proportional to c m Cd^Kir/Mg, which 
is generated through 778—771 mixing. 
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together with the relations ( 2.12|) , immediately implies Lf = 2Lf, independent of the 



number of scalar resonances. Within the uncertainties, this relation is fulfilled by the most 
recent determination of the chiral couplings L5 and L§ |43|| , 

L 5 = (0.91 ± 0.15) ■ 10~ 3 and L 8 = (0.62 ± 0.20) ■ 10~ 3 , (2.13) 

suggesting that both the resonance saturation as well as the short-distance constraints are 
reasonable approximations. However, we only have three linearly independent relations for 
the four scalar couplings q, c m , c' d and c' m , so that we have to make further assumptions 
to obtain estimates for these couplings. 

The estimate ( |2.3| ) with only one scalar resonance led to q = c m and thus it seems 
plausible to keep this constraint in addition. One immediate consequence of the constraint 



and the second relation (|2lf ) is c' d = c' m . Employing the first of the short-distance con- 
straints (|2^ ) and the relations fl2.12|) , the remaining couplings c d and c' d can be calculated 
with the result: 

Cd = Cm = 37.0 MeV and c' d = c' m = 27.7 MeV. (2.14) 

Here, the phenomenological value (|2.13|) for L 5 , as well as Ms = 1.4 GeV and M' s = 
1.9 GeV have been used. The approach of utilising the short-distance constraints, together 
with Cd = c m has also been followed in our fit (6.10) of ref. |l|. The fit then resulted in 

Cd = Cm = 23.8 MeV and c' d = c' m = 39.6 MeV. (2.15) 

These values would correspond to Lf = 0.72 ■ 10 -3 , which in view of the phenomenological 



result (|2.13|) appears a little low, although nevertheless acceptable. One should note that 
at the tree-level there is some ambiguity in the scalar resonance mass Mg, and e.g. taking 
M s = 1.2 GeV, eq. ( gTg) would change to c d = 27.6 MeV and c' d = 37.0 MeV, much closer 



to the result ( p,15|) . 



A different, although also not implausible assumption would be universality of the chiral 
couplings, namely c d = c' d and c m = c' m . These constraints allow us to calculate Lf from 
the first of the short-distance relations: 



= + = ' 84 ' 10 "' (2 ' 16) 

surprisingly consistent with the phenomenological value (|2.13|) . On the other hand, just 
using resonance saturation of the chiral couplings Q2.12|) and their values (|2.13|) , yields 

Cd = c ' d = 29.1 MeV and c m = c' m = 39.7 MeV, (2.17) 
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which also satisfies the first short-distance relation rather well. However, these parameters 
violate the second of the short-distance constraints. Enforcing the second constraint in 
addition again requires Cd = c m , and we fall back to the first assumption. The values of 
the chiral scalar couplings which correspond to this last scenario are 

c d = c m = c' d = c' m = 32.7 MeV. (2.18) 

The ranges of the scalar couplings which arose in the discussion above should give some 
indication of their present uncertainties which we conclude to be about 30%. 

3 Analyticity and unitarity constraints 

In this section, we present the formalism employed to calculate the coupled channel Ktt, 
Kr] and Krj' scalar form factors up to around 2 GeV. The results (|2.7| ) and ( |2.8| ) are 
only valid at low momentum transfers. We have been able to resum the local leading (in 
1/iVc) contributions to all chiral orders, through the resonance propagators. Analyticity 
and unitarity allow us to obtain further constraints on the scalar form factors. We shall 
first derive unitarity relations obeyed by the scalar form factors which link these quantities 
to the S-wave 1=1/2 Ktt, Krj and Krj' partial wave amplitudes. These amplitudes will 
be taken from our previous work where they were studied in detail up to 2 GeV in the 
framework of xPT with resonances, supplemented with a suitable unitarisation procedure. 

3.1 Unitarity relations 

Unitarity of the scattering matrix S immediately implies the following identity for the 
T-operator, 

T-T f = iT-T\ (3.1) 

where T is defined by S = 1 + i T. To obtain the form factors in question, the previous 
general relation is applied to the transition of the states \K^>k) (where 0i = tv, <p2 — V and 
03 = rj') to the vacuum state. The scalar operator giving rise to this transition should have 
strangeness \S\ = 1 and isospin / = 1/2. Let us remark that we are considering a pure 
strong interaction problem without electroweak corrections. 

Inserting a complete set of intermediate states on the right-hand side of eq. fl3.1|) , and 
restricting this set to the two-particle states \K(pi), we arrive at 

(0\T\K<l> k ) - <O|Tt|#0 fc ) = iJ2 6 ( a (0|T|^) [(K^lK^dcose, 

i cmy s J 

(3.2) 
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where s = p\ with pi being the total four-momentum of the state \K<f>i), qi = \k</>J(2\/s) 
is the modulus of the centre-of-mass three-momentum for this state, and 9 is the corre- 
sponding scattering angle. The function Xtcfa is defined in appendix A, and the trivial 
centre-of-mass motion has been removed. 

Performing a partial wave decomposition of T tk (s,cos8) = (K<j)i\T\K(j)k), we have: 

oo 

T k (s,z) = 167r£(2J + l)fj*(s)^(z), (3.3) 
1=0 

where Vi(z) is the Legendre polynomial of degree I with I being the orbital angular mo- 
mentum and t\ k (s) are the partial wave amplitudes. Making use of time reversal invariance 
and inserting the decomposition ( |3.3j ) into eq. (|3~2|), only the S-wave component survives 
the angular integration. We thus obtain our central unitarity relation for the scalar form 
factors F k {s) = (0\T\K</> k ): 

lmF k (s) = J2^)Fi{s)ti k (s)*, (3.4) 

i 

where, for convenience, we have defined the quantity a*(s) = 9(s — s t hi) 2qi(s)/yfs. For an 
appropriate choice of the scalar T-operator with isospin 1/2, in accordance with eq. Ql.lj) , 
the form factors F k (s) are related to fx(s) of section 1 by: 

C C i 

Fktt(s) = ficM, F Kr) (s) = -jSf Kv (s) } F Kv ,(s) = -^-f Kr)/ (s). (3.5) 

We use an isospin basis of states \K(f>k). The coupling of the \Ktt) state with I — I 3 — 1/2 
to the scalar current in (1.1) is a factor v3 larger than the one of |i^ + 7r°). Too keep the 
same normalisation as for /_K+ 7r o(s), we rescale all form factors by a global factor l/y/3. 

As stated above, the sum in eq. ( |3.4| ) is restricted to two-particle intermediate states. 
Thus some comments about the role of multiparticle states are in order. The lightest 
multiparticle state contributing to this sum is \K7r7nr). From the theoretical side, its 
contributions are suppressed both in the chiral and large-iV c expansions. Experimentally, 
in [j44|, [45H it has been established that the 1=1/2 S-wave Kir amplitude is elastic below 
roughly 1.3 GeV with Krj' being the first relevant inelastic channel. These conclusions have 
been confirmed in |1] where it was demonstrated that the contribution of the Kr] channel to 
the S-wave amplitude could be neglected. However, it was also found in this work that for 
energies higher than about 2 GeV other inelastic channels become increasingly important. 

The unitarity relation ( |3.4| ) given above poses tight constraints on the scalar form 
factors. For the elastic case one simply has: 

ImFi(s) = a 1 {s)F 1 {s)tl 1 ( S y (3.6) 
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which implies the well known Watson final state theorem stating that the phase of 
Fi(s) is the same as the one of ^(s). This is obvious from the previous equation since its 
left-hand side is real. 

Let us now discuss the two-channel case, which is very appropriate since the contribu- 
tions from the Kr\ state are negligible to a very good approximation. Substituting ImFj.(s) 
by (-Ffc(s) — F k (s)*)/2i in eq.( |3.4D , one finds above the threshold of the Krf: 

F x {s) = (l + 2ia 1 (s)tl 1 (s))F 1 ( S y + 2t a 3 (s)tl 3 (s) F 3 (s)* } (3.7) 

F 3 (s) = 21^(3)^(3) F 1 (s)*+(l + 2ia 3 (s)tf( S j)F 3 (sy, (3.8) 

where, although working with two channels, we have used the label 3 to indicate the Krf 



state as introduced before. Following ref. |47[], we now express F 3 (s) in terms of -Fi(s) from 
eqs. (|3) and Q. 

The partial wave amplitudes t™ n can be parametrised by introducing the symmetric 
and unitary 2x2 S-matrix S mn = 5 mn + 2i^fa m a n t™": 



S 



Tj exp 2iSi(s) 
i^l - r\ 2 exp i(5 1 (s) + 5 3 (s)) 



- ?7 2 expi(5i(s) + S 3 (s)) 
r] exp 2iS 3 (s) 



(3.9) 



with 7] = cos 2a (sin 2a = ^1 — r/ 2 ) the inelasticity parameter (0 < r] < 1) and Sk(s) the 
phase shift of channel k. We further express the form factors as Fk(s) = fkexpi(Sk + 4>k) 
with fk and (fik real functions and fk > 0. Taking real and imaginary parts in eqs. ( |3.7| ) 
and (|3.8| ), the following set of relations emerges: 



[1 — cos 2a) cos0i/i 
(1 + cos 2a) sin 4>xf\ 



/03 



sin 2a sin0 3 / 3 , 
sin 2a cos0 3 / 3 ; 



(3.10) 
(3.11) 



plus two analogous relations with the labels 1 and 3 exchanged. The latter relations are, 
however, linearly dependent to the first two. Dividing the two eqs. (|3.10|) and (|3.11|) results 



in: 



tan <pi tan 3 
whereas adding them in quadrature yields: 



tan a , 



tan 2 a + (cot 2 a — tan 2 a) sin 2 <pi . 



(3.12) 



(3.13) 



Thus, once fi and <pi are known, / 3 and <p 3 can be calculated from eqs. ( |3.12|) and (|3.13| ) 
The treatment in the full three-channel case will be discussed further below. 
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3.2 Dispersion relations 

The scalar form factors Fk(s) are analytic functions in the complex s-plane, except for a 
cut along the positive real axis, starting at the first physical threshold s t hi = (Mk + M^) 2 , 
where their imaginary parts develop discontinuities. They are real for s < s t hi- As should 
be clear from eq. ( |3.4|) , their imaginary parts just correspond to the contributions from all 
possible on-shell intermediate states.0 

Analyticity relates the real and imaginary parts of Fk(s) through dispersion relations: 



1 °r ImF ( s') 

F k (s) = - Ids'— kK '— + subtractions . (3.14) 

TT J (,$' — s — iO) 



Hkl 



From unitarity, as shown above, ImFk(s) takes the form of eq. (|3.4j ). With the reasonable 
assumption that for s going to infinity, Fk(s) vanishes sufficiently fast, the form factors 
satisfy dispersion relations without subtractions. This assumption was already explored in 
section 2. 

Because Kir scattering is basically elastic up to around 1.3 GeV, it is again instructive 
to consider the single channel case. Then the dispersion relation takes the form: 



flW . I L«M. (3 . 15) 

K ' if J (s'-s-iO) K ' 

Hh l 

In this case the partial wave amplitude can be expressed as cr^ 1 = sin 5i exp(i5i), and the 
previous eq. ( |3.15| ) admits the well known Omnes-exponential solution [fi8|| : 



F 1 (s) = P(s) exp - / da' „ , U ; - , (3.16) 




with P(s) being a real polynomial to take care of the zeros of Fi(s) for finite s and 5i(s) 
is the S-wave I = 1/2 elastic Ku phase shift. 

In writing eq. ( |3.16| ), we have included one subtraction at the origin since generally 
S\(s) tends to a constant for s going to infinity. A general solution, valid for any number n 
of subtractions at an arbitrary subtraction point Sq, has been given in refs. where 



the equivalence of the different expressions has been discussed in detail. The form factor 
F\(s) in the elastic case will be discussed further in our numerical analysis below. 

In order to implement the dispersion relation for coupled channels, we neglect, for the 
moment, the Krj contribution and only consider the more important Kit and Krj' channels. 
2 We are excluding the presence of bound state poles below the threshold of the Kir state. 
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For this case, two coupled integral equations arise: 

v ; 7T 7 (s'-s-iO) Ti J (s'-s-iO) 

s thl s th3 

= 1 L«WW t l L*«M. (3,7) 
v; vr 7 (s'-s-iO) vr J (s'-s-iO) v 7 

The integral equations ( |3.17|) will be solved iteratively according to a procedure already 
applied in ref. [50]. In the first step initial functions Fjf\s) for the form factors, to be 



specified further below, are inserted in the right-hand side of eqs. ( |3.17| ), and new resulting 
form factors Fjp(s) are calculated. Then the procedure is iterated until after n steps the 
procedure converges and the final scalar form factors Fj, n \s) are obtained. 

The dispersion relations for the three-channel case including the Krj contribution are 
completely analogous to the eqs. ( |3.17| ). Thus the explicit expressions have been relegated 
to appendix B. As for the two-channel case, they are solved iteratively. Further details on 
the numerical methods will be given in the next section below. 



4 Numerical analysis 

In this section we present the different solutions to the dispersion relations of eqs. (|3.15| ), 
( p,17| ) and ( P-l| ) for the scalar Kn, Krj and Krf form factors. The required partial wave 



amplitudes are taken from our previous work ||. First, we investigate the single channel 
elastic case for which the analytical Omnes solution is available. Then the inelastic two- 
and three-channel problems are treated for which we have to resort to numerical methods. 

4.1 Elastic Kit channel 

In the elastic case, the analytical Omnes solution is given by eq. (|3.16|) . It only depends 
on the elastic Kit phase shift ^i(s). However, we still have to fix the polynomial ambiguity 
P(s). This can be achieved from the assumption that Fi(s) should vanish for s going to 
infinity. 

Investigating Fi(s) in this limit, from eq. fl3.16Q one finds: 



lim Fi(s) = e i5l °° lim P(s) ( — ) , (4.1) 

s— >oo s^oo \ S J 

with <5 loo = lim^oo <5i(s). In the single channel case with one resonance, 5^ is expected 
to be equal to n [^TJ. Assuming Fi(s) to vanish at infinity thus requires that P(s) should 
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be constant. This constant can be determined from the normalisation of i*i(s) at s 
Then the Omnes formula takes the form: 



0. 



F 1 (s) = F 1 {0)exp[- f da' - W) 

\ 7T J S'[S' — S — 



Sthl 



zO) 



(4.2) 



In our numerical analysis, for -Fi(O) we will use the next-to-leading order xPT result of 
eq. (|2.7f) , namely Fi FT (0) = /^(O) = 0.981. This is very close to the leading order value 
1, demonstrating that at s = 0, as expected, xPT works extremely well. 
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Figure 1: |-Fi(s)| in the elastic case. Long-dashed line: xPT fit of eq. (4.9) of solid 
line: ;\PT fit of eq. (4.10) of |Q; short-dashed line: K- matrix fit of eq. (4.13) of dotted 
line: next-to-leading order xPT result of eq. 



Inputs for the phase shift Si(s), required in eq. ([4.2|) , will be taken from different fits 
of our previous work on S-wave Kit scattering |l|. In fig. [1] the absolute value of the form 
factor, |Fi(s)|, is displayed for three different elastic fits, as a function of the invariant mass 
of the Kit system M KlT = yfs. The long-dashed line corresponds to the xPT fit of eq. (4.9) 
of HJ, the solid line to the fit of eq. (4.10) of |l[], and finally the short-dashed line 
corresponds to the K-matrix fit of eq. (4.13) of M. For comparison, as the dotted line we 
also show the next-to-leading order %PT result of eq. ( |2.7| ). This comparison shows that 
up to the Ktc threshold region, next-to-leading order ;\TT gives a good description of the 
scalar form factor. 

As is obvious from fig. [I], although all three fits give very good representations of the 
Kit scattering data in the elastic region, the corresponding form factors, especially in the 
resonance region, behave very differently. Let us discuss these findings in more detail. In 
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the resonance region the form factor corresponding to our best elastic fit of eq. (4.10) of Jl| 
is almost flat, displaying no significant resonance structure. This can be traced back to 
the fact that for large s, 5j 4 ' 10 ^(s) tends to zero, whereas for the passage of one resonance 
the phase shift should go to n. The fit of eq. (4.9) of |l|], including the short-distance 
constraint discussed in section 2, on the other hand, displays this behaviour. However, in 
this case the scalar form factor in the resonance region is rather large. This is due to the 
fact that for this fit the width of the ifg(1430) comes out to be very small, about a factor 
of three below the experimental average. (See section 7 of 0.) Nevertheless, one has to 
strongly stress that a proper study of the resonance region can only be performed after 
including the Krj' with a threshold around 1.45 GeV, as shown below. 

With respect to the general behaviour discussed above, the form factor corresponding 
to the K-matrix fit of eq. (4.13) of [|I|] appears to be the most realistic one. 5^' 13 \s) tends 
to 7i for large s and also the resonance parameters of the -fTg(1430) for this fit are close 
to the experimental average. The K-matrix fit is also different from the chiral fits because 
the corresponding ansatz for the partial wave amplitude is free of left-hand cuts. In this 



case it is even possible to give a closed expression for the form factor [47. 51 



F,ls) - n (1 " s/Sw) Fl(0) U3) 

— 1/2 

where C(s) = 16tcJkit(s) and the explicit expression for the K-matrix ansatz K (s) 
is given in eq. (4.12) of [![]. The s Pi and s Zi are the locations of the poles and zeros 
of (1 — C(s)K(s))^ 1 which have to be removed in the form factor. In our particular 
case we have one pole at s Pl = —2.011 GeV 2 and two zeros at s Zl = —5.821 GeV 2 and 
s Z2 = It is an easy matter to verify that the representation ( |4.3|) yields the same 

scalar form factor as the Omnes formula 01. 2D . 



4.2 The two— channel case 

We now turn to the two-channel case, for which the important effects of the Krj' are 
considered in the coupled channel integral equations ( |3.17| ). Let us first briefly describe 
our numerical treatment of these integral equations. 

The first step consists in rewriting the coupled integral equations according to the 
procedure described in appendix C.l. As central values for the two parameters which are 
introduced in this context we take s cu t = 9 GeV 2 and b = — 0.8. We have, however, varied 
these parameters to test the numerical stability of our approach.f] Next, the resulting 
3 Onc of us, J.A.O., has implemented a different numerical approach, finding agreement for the resulting 
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integral equations are solved iteratively according to a procedure already employed in 



ref. [50 



For the iterative solution we require initial functions Fi°\s) and F^°\s). For 



we use the form factor which results from employing the Omnes formula (f4.2|) . The phase- 
shift 5i(s) is taken according to the decomposition ( |3.9|) . As a starting point, the form 
factor i*3 (s) can be set to zero. We have verified that other initial choices for this function 
lead to the same final form factors. 

Since the dispersion integrals are calculated according to the principal value prescrip- 
tion, the immediate results of our first iteration step are the real parts ReF^ (s) and 
ReF^^s). From these results, together with the initial imaginary parts ImFf and ImFf, 
and making use of the central unitarity relation (|3.4Q , the new full form factors and F^ 
are obtained. This procedure is repeated, until after n steps it has converged reasonably 
well and the final form factors F^ and F 3 can be extracted. 

In what follows, we shall explicitly investigate the form factors corresponding to the 
two-channel fits (6.10) and (6.11) of ref. [[jj to the S-wave f^7r-scattering data. In a 
forthcoming publication [[5^], these form factors and the corresponding scalar spectral 
function will be utilised to calculate the mass of the strange quark in the framework of 
QCD sum rules. 

The absolute value of the resulting inelastic two-channel Ktt and Krj' form factors 
Fi(s) and F 3 (s) is displayed in figures || and |3|. The solid line corresponds to the x?T fit 
of eq. (6.10) of ref. M and the long-dashed line to the fit (6.11). Depending on the high- 
energy behaviour of the phase shifts Si(s) and Ss(s), the dispersion relations admit one 
or two linearly independent solutions [j47H . For the fits (6.10) and (6.11), with s going to 
infinity, Si(s) tends to tt whereas Sz{s) tends to zero. In this case there is only one solution 
to the two-channel dispersion relation. However, we still have the freedom to normalise 
both solutions Fl(s) and F 3 (s) by a common factor. We have used this freedom to again 
fix F* PT (0) = f Kn (0) = 0.981. 

The value of the Kr( form factor at zero momentum is then determined and we obtain 
F 3 6 ' 10 ^(0) = 0.430 and F 3 ^(O) = 0.409. These results should be compared with the tree- 
level expectation from U(3)l x U(3)# in the large- N c limit including resonances of 
eq. (|J), F 3 XPT (0) = 2^2/3/^(0) = 0.696. We observe that the resulting value for F 3 (0) 
from the solution of the coupled channel dispersion relation turns out to be somewhat 
lower than the tree-level result. We shall come back to a discussion of this point below. 

For comparison, as the dotted line in figure ^|, we have displayed the form factor that 
form factors. 
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Figure 2: |i*i(s)| in the two-channel case. The solid and long-dashed lines correspond to 
the xPT fits of eqs. (6.10) and (6.11) of [l] respectively. For comparison, the elastic K- 
matrix fit of eq. (4.13) is displayed as the short-dashed line, and the dotted line represents 
the Omnes formula (|4.2|), evaluated with S^' 10 \s). 



results from evaluating the Omnes formula eq. (|4.2[) with 5± 



,(6.10), 



[s] has been 



used as the starting point for the iteration of the integral equations. From figure it is 
clear that, apart from the region of the second resonance, the Omnes representation of 
Fi(s) already gives a reasonable description of this form factor. In addition, as the short- 
dashed line in figure we have also plotted |.Fi(s)| corresponding to the elastic K-matrix 
fit of eq. (4.13) of (![]. Here we observe that the corresponding form factor behaves rather 
differently in the Kir threshold region. This should come as no surprise, because as has 
been already discussed at the end of section 4 of ref. [0], also the scattering lengths for this 
fit come out very different compared to x?T. 

As was already mentioned above, for the fits (6.10) and (6.11) the total phase motion 
of Si(s) + £3(5) for s going to infinity only reaches 7r, whereas on general grounds for two 
resonances it is expected to approach 2ir. The reason for this behaviour of our unitarised 
chiral fits can be traced back to a deficient description of the experimental data [4~5| above 
roughly 1.9 GeV. To improve the description in the region of the second resonance and 
above, we have performed new fits with a K-matrix ansatz. In the fitting process, the 
K-matrix ansatz is matched smoothly to the chiral fits at an energy around 1.75 GeV. 



15 



5.0 
4.0 





II 


I 


1 1 1 1 1 1 














/ v 




■ 








1 


\\ 










1 


\\ 










1 1 


\\ 










r 


\\ 










i 


\\ 










/A 


\\ 










/ 


\\ 
\\ 
\ 






I 












1.0 
0.0 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

M Kn [GeV] 

Figure 3: |i*3(s)| in the two-channel case. The solid and long-dashed lines correspond to 
the xPT fits of eqs. (6.10) and (6.11) of JI| respectively. 

For the K-matrix we use a resonance plus background parametrisation: 

lA ) ~ (MJ, - s) + [1 + (s/cY) ' {AA) 

with k = 0, 1 if bij = and k = 2 otherwise. Except for the first case k = 0, this K-matrix 
is constructed such that it vanishes linearly at large s, and above 1.75 GeV it replaces 
the matrix N(s) in eq. (6.1) of ref. IJ. Because of time-reversal invariance, K(s) has to 
be symmetric. Above 1.75 GeV, we then fit this ansatz to the data set A of |15] which 
extends up to 2.52 GeV. Simultaneously a smooth matching to either one of the unitarised 
chiral fits (6.10) or (6.11) of |j] at around 1.75 GeV is imposed in the fitting procedure as 
well. To judge the dependence of the K-matrix fits on the background parametrisation, we 
have actually calculated four different types of fits. All fits have been performed with the 
program Minuit p3| , and the parameters of these fits are compiled in table |I|. As examples, 



in figure [| we display the experimental Ktt scattering data of refs. [44,^5|, together with 
the unitarised chiral fit (6.10) of M (solid line) and the new K-matrix fits with improved 
behaviour in the region of the second resonance (6.10K1) (dashed-dotted line) and (6.10K4) 
(dotted line). The corresponding curves for the other fits look very similar. 

For the first type of fits (6.10K1) and (6.11K1), as the background we have just taken 
constants a^. This entails that the T-matrix at large s does only fall off as l/ln(s), and 
that the total phase motion at infinity, 5 loo + <5 3oo is only tt, analogous to the case of the 
unitarised chiral fits. For the convenience of the reader, the separate phases at infinity, 5^ 
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Fit 


(6.10K1) 


(6.10K2) 


(6.10K3) 


(6.10K4) 


(6.11K1) 


(6.11K2) 


(6.11K3) 


(6.11K4) 


Mo [GeVl 

a L J 


2.040 


2.040 


1.904 


1.909 


1.837 


1.837 


1.838 


1.837 


r x [GeV] 


1.007 


1.003 


0.567 


0.617 


0.495 


0.501 


0.449 


0.433 


ro [GeVl 


1.420 


1.406 


0.550 


0.648 


0.595 


0.603 


0.499 


0.474 




-0.723 


-0.801 


-0.703 


0.221 


-0.350 


-0.407 


-0.647 


-0.346 




-0.905 


-0.990 


0.495 


1.797 


-0.356 


-0.415 


0.565 


1.741 


^22 


-1.450 


-1.580 


-1.276 


-0.159 


-0.567 


-0.656 


-0.390 


0.750 


6n [GcV- 2 ] 








0.106 


-0.308 








0.145 


0.028 


6i a [GcV- 2 ] 








-0.214 


-0.738 








-0.224 


-0.606 


fe 22 [GcV- 2 ] 








0.366 


-0.123 








0.055 


-0.326 


c [GeV 2 ] 


oo 


25 


25 


4 


oo 


25 


25 


4 


X 2 /33 dof 


74.9 


75.5 


42.3 


45.2 


83.3 


87.4 


43.2 


38.0 


^loo; 03oo 


7V, 


71", 7T 


2vr, 


2tt, 


7V, 


TV, 71 


71, TX 


2tt, 



Table 1: Different K-matrix fits corresponding to the fits (6.10) and (6.11) of JT| in the 
region above 1.75 GeV. For a detailed discussion see the text. 

and <5 3oo , have also been collected in table [L| For these two fits, again F 3 (0) is determined, 
and comes out as F 3 (6 - 10X1) (0) = 0.591 and F 3 (6 - 11/o) (0) = 0.574 for (6.10K1) and (6.11K1) 
respectively, much closer to the tree-level result F* PT (0) = 0.696. As was already remarked 
above, at the tree-level the scalar resonance mass Ms has some ambiguity, and e.g. using 
Ms = 1.2 GeV in eq. fl2.5|) would yield F 3 XPT (0) = 0.586, in perfect agreement with the fit 
results for -F 3 (0). 

To further investigate the quality of our fits, in the following we shall include additional 
low-energy constraints on the scalar Ktt form factor which come from pure next-to-leading 
order %PT Rather precise information on Fi(s) is available at the so called Callan- 
Treiman point M\ — M%, because at that point Fi(A Kn ) only differs from F K j 'F n by terms 
of order m u or m d , 

Fi(Akit) = ^ + A CT , (4.5) 

where numerically Act = — 3.0 • 10~ 3 at the next-to-leading order in xPT, and the analytic 
expression can be found in the second of refs. 0]. Taking most recent values for the decay 
constants F% and F n from ||54fl , we conclude that ^(A^^) = 1.22 ± 0.01. 



Furthermore, from x?T we also have information on the slope of the Ktt form factor 
at zero energy, Di(0) = F 1 '(0)/F 1 (0). In at the next-to- leading order, one obtains 
-Di(O) = 0.873 GeV -2 , where the uncertainty due to higher order corrections was estimated 
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Figure 4: Two-channel fits for ao and 0o- The experimental data are given by: solution B 
of [44 1 full circles; solution A of [45| full squares. The solid line represents our unitarised 
chiral fit (6.10) of ref. [HJ, and, as examples, the dashed-dotted and dotted lines show the 
new K-matrix fits (6.10K1) and (6.10K4) respectively, which display an improved behaviour 
in the region of the second resonance. 



to be around 20% 0. This value can also be compared to experimental data. The slope 
of the scalar form factor has been measured in charged as well as neutral Kaon decays, 
with the results A = 0.006 ± 0.007 and A = 0.025 ± 0.006 respectively |54|. Here, the 
parameter Ao is defined by Ao = -Di(0)M 2 . Taking the average of the experimental results 
and enlarging the error according to the standard procedure of the Particle Data Group 
for inconsistent measurements, we obtain Ao = 0.017 ±0.009, whereas from the xPT result 
together with the isospin average M n = 138 MeV, we also find A = 0.017, displaying 
complete agreement for the central results. 

Let us now compare -Di(O) and especially Fi(Aktt) with the numbers emerging from 
our fits. For the K-matrix fits which accomplish a good description of the high energy 
data above 1.75 GeV, we get F-f 6 ' 10 ^ ( A K7r ) = 1.220 and F^' 11K1 \A Kn ) = 1.219, in 
impressive agreement with the precise %PT expectation of eq. ( }4.5| ). In addition, the slope 
is found to be L>f 10in) (0) = 0.894 GeV" 2 as well as L>f 11K1) (0) = 0.893 GeV" 2 , also 



very close to the ;\PT value. On the other hand, for the original unitarised chiral fits 

(6-10)/ A .. ^ = 1.262 and (A K7T ) = 1.260, as well as 

0.909 GeV" 2 . Although the values for £>i(0) 



the corresponding results are F± ' (A^ 



Dj 6 - 10 ^) = 0.918 GeV"" and Df^O) 



,(6.11), 
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are compatible with our previous results, Fi(Akk) turns out too large for both fits, again 
reflecting an insufficient description of the higher energy region for these fits. 



Fi(A^) 


F 8 (0) 


Di(0) [GcV- 2 ] 




F 3 (0) 


L>i(0) [GcV" 2 ] 




(6.10K2) 








(6.11K2) 




1.21 


0.636 


0. 


348 


1.21 


0.612 


0.850 


1.22 


0.587 


0. 


363 


1.22 


0.565 


0.866 


1.23 


0.538 


0. 


379 


1.23 


0.518 


0.882 




(6.10K3) 








(6.11K3) 




1.21 


0.586 


0. 


328 


1.21 


0.618 


0.777 


1.22 


0.567 


0. 


351 


1.22 


0.587 


0.805 


1.23 


0.547 


0. 


374 


1.23 


0.556 


0.834 




(6.10K4) 








(6.11K4) 




1.21 


0.622 


0. 


351 


1.21 


0.615 


0.858 


1.22 


0.590 


0. 


368 


1.22 


0.579 


0.874 


1.23 


0.557 


0. 


386 


1.23 


0.544 


0.890 



Table 2: -F 3 (0) and -Di(O) for the unitarised chiral plus K-matrix fits (6.10K2-4) and 
(6.11K2-4) with 8 loo + 8 3o o = 2tt, chosen such that Fi(A x?r ) = 1.22 ± 0.01. 

For the previous fits, the total phase motion 5 loo + 8 3oo was only n, although for two 
resonances 2n would be expected. The differences are due to the background accompanying 
the bare K-matrix poles in eq. ( |4.4| ). The latter behaviour can be achieved by forcing the 
K-matrix to vanish at infinity as will be implemented in all our remaining fits. In this 
way, we also test the stability of our results under changes in the parameterisations that 
induce different T-matrices well above the second resonance region around 2 GeV. As a 
first step, in the fits (6.10K2) and (6.11K2), we introduce an additional cutoff parameter 
c = 25 GeV 2 . As long as c is large enough, the sensitivity on this parameter is rather 
small and it cannot be incorporated as a fit parameter. As can be seen from table [TJ, 
the remaining fit parameters do depend only little on this modification. Nevertheless, 
now both phase shifts S\ and 8 3 approach it at infinity. As a consequence, the dispersion 
relations ( ^.17| ) admit two linearly independent solutions. Thus, besides F 1 (0), now we are 
in a position to also fix 1*3(0) as an initial condition. In our opinion the best information 
available is Fi(Aktt), and therefore we fix Fs(Q) such as to fullfill the constraint on this 
quantity presented above. 

Our results for this exercise are shown in table Before we come to a detailed discussion 
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of these results, however, let us first present the remaining fits of table ffl. As can be read 
off this table, the \ 2 of the K-matrix fits (6.10K1), (6.10K2), (6.11K1) and (6.11K2) is 
larger than 2 for all cases. Let us remark that this x 2 only corresponds to the data points 



of [45] above 1.85 GeV, which have not yet been included in the fits (6.10) and (6.11), and 
the matching condition around 1.75 GeV. To improve the quality of our K-matrix fits, 
we have decided to include an additional term in the background. To maintain the 
required falloff of the T-matrix, for these fits the parameter k = 2. In order to study the 
cutoff dependence we have performed this type of fits with the two values c = 25 GeV 2 
and a lower cutoff c = 4 GeV 2 . 

For the fits (6.10K3), (6.10K4) and (6.11K4), we observe that the phase shifts at infinity 
turn out to be = 2n and S^oo — 0. The reason for this change in the behaviour of the 
phase shifts lies in the fact that the second resonance couples almost equally to the first and 
the third channel, so that a small change in the fit parameters can shift the phase motion 
of 7r from the third to the first channel. For the parameter set where this happens, we also 
observe a vanishing inelasticity and thus zeros in the diagonal elements of the S-matrix. 
Therefore the phases 5i(s) and ^(s) at this point become ambiguous, although their sum 
must be continuous (modulo 2ir) since the matrix element S\% of eq. (|3.9|) is non-vanishing. 
As far as the phase shifts at infinity are concerned, the fit (6.11K3) behaves as (6.10K2) 
and (6.11K2). However, this fit displays an artificial rapid phase motion at an energy 
around 13 GeV, and although it has a rather good \ 2 an d the results for F 3 (0) in table || 
turn out reasonable, we shall discard it in the following. 

Our results for the absolute values of the form factors F\{s) and F$(s) in the case of the 
unitarised chiral plus K-matrix fits are shown in figures [5] and |[ The dotted lines represent 
the fits (6.10K1) and (6.11K1), the solid lines (6.10K2), (6.10K3) and (6.10K4), and finally, 
the dashed lines correspond to (6.11K2) and (6.11K4). Since all form factors for these fits 
are rather similar we have not tried to discriminate the curves further. Nevertheless, for 
|F 3 (s)|, like in figure |3|, there are differences in the form factor between the fit (6.10) and 
(6.11) which are most prominent around 1.2 GeV. These differences can already be seen 
in the T-matrix and they tend us to consider the fit (6.10) more physical, as was already 
concluded in ref. ffl. 

We shall return now to a discussion of table |. In figures [7| and we display | _F\ (s) | 
and 1^3(5)1 for our fits (6.10K3) and (6.10K4), while varying ^(A^). The dotted lines 
correspond to Fi(A^) = 1.21, the solid lines to 1.22, and the dashed lines to 1.23. The 
curves which are higher at the _K"q(1430) resonance always represent the fit (6.10K4). As 
can be seen from table §, for smaller values of Fi(Aktt), Fs(0) comes out larger and vice 
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Figure 5: in the two-channel case with K-matrix ansatz. The dotted lines corre- 

spond to the fits (6.10K1) and (6.11K1), the solid lines to (6.10K2), (6.10K3) and (6.10K4), 
and the dashed lines to (6.11K2) and (6.11K4). 

versa. Averaging our results for the K-matrix fits, for F 3 (0) we then extract 

F Krj/ (0) = 0.58 ±0.04. (4.6) 

As the uncertainty we have chosen the maximal variation for the fits K3 and K4 with the 
better x 2 - The slightly larger variations in the case of K2, to our minds are due to the 
worse description of the experimental data. Also the derivative of the form factor i*i(s) at 
zero, -Di(O), comes out completely consistent with our expectation from xPT, with some 
variation for the fit (6.11K3), which anyway has been discarded for other reasons, although 
even this variation stays within the uncertainties. Finally, it is worth to note here that the 
predicted values for F 3 (0) from the fits (6.10K1) and (6.1 1K1) are in perfect agreement 
with the previous central value of F 3 (0), and their difference is completely accounted for 
by the estimated error in eq. (|4.6|) . 

4.3 The three— channel case 

In the following, we investigate the solutions of the dispersion relations ( P-l|) in the full 
three-channel case for our unitarised chiral fits (6.10) and (6.11) of ref. JT|. In the three- 
channel case, we have not improved our fits in the region of the second resonance, because 
on the one hand just from the Kir scattering data there is not enough information to fit all 
parameters of a three-channel K-matrix, and, on the other hand, as was already mentioned 
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Figure 6: | -F3 ( s ) | in the two-channel case with K-matrix ansatz. The dotted lines corre- 
spond to the fits (6.10K1) and (6.11K1), the solid lines to (6.10K2), (6.10K3) and (6.10K4), 
and the dashed lines to (6.11K2) and (6.11K4). 

above, the effects of the Krj channel will turn out so small that the more elaborate two- 
channel analysis with only Ktt and Krj' for our purposes is completely sufficient .0 As we 
shall verify below, this also entails that the form factors Fi(s) and F 3 (s) are completely 
stable under the inclusion of the scalar Krj form factor Fz(s). Therefore, this form factor 
will play a negligible role in the determination of the mass of the strange quark |52|] . 

The numerical solution of the dispersion relations proceeds along the same lines as was 
already discussed in the last section for the two-channel case. Again, as initial functions 
for the form factors we can take the Omnes solution for F±°\s), whereas the form factors 
F2°\s) and F 3 \s) can be set to zero. Nevertheless, in order to prove the stability of the 
two-channel solutions (6.10) and (6.11) under the inclusion of the Krj channel, we have also 
used these solutions as initial functions Fi°\s) and F 3 \s) in the present three-channel 
case. Then the solution of the dispersion relations are iterated, until the procedure has 
converged and the form factors F 2 (s) and F 3 (s) can be extracted. Of course, with 

both choices, the same final form factors are obtained. 

In figure |9|, we display the absolute values of the form factor F±(s) for the three-channel 

4 We have not taken into account the fit (6.7) of ref. 0, because when the Krj channel is removed for 
this fit, the x 2 becomes rather large, increasing by almost a factor of two. On the theoretical side, this fit 
does not fullfill the short distance constraints presented in section 2. The strong sensitivity of the fit (6.7) 
on the Krj channel, and also the very small value of c m , tend us to disregard this fit for further analysis. 
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Figure 7: |Fi(s)| for the fits (6.10K3) and (6.10K4) while varying Fi(A Kv ). The dotted 
lines correspond to Fi(Ax n ) — 1-21, the solid lines to 1.22, and the dashed lines to 1.23. 
The curves which are higher at the ^(1430) resonance always correspond to (6.10K4). 



fits (6.10) and (6.11) of ref. ||], and in figure [K] the corresponding form factors |-p2(s)| and 
|i 7 3(s) | are shown. Like in section 4.2, the notation of the lines is (6.10) solid line and 
(6.11) long-dashed line. The three-channel form factors are plotted as thick lines and for 
comparison, in addition, as the thin lines we have displayed the two-channel form factors 
for the fits (6.10) and (6.11). In figure 10, the two lower lines correspond to |F 2 (s)|. As can 
be observed from this figure, the form factors F 2 (s) turns out much smaller than F 3 (s), and 



for the scalar spectral function, further discussed in ||52|| , it will be completely negligible. 
Furthermore, the two- and three-channel form factors are found to be very similar which 
again supports the finding that the Kr) channel is rather unimportant. 

Like in the case of the fits (6.10) and (6.11) for two channels, there is only one solution 
to the dispersion relations and we just have the freedom to fix the normalisation i*i(0). 
The form factors ^(0) and ^(0) are then determined from the dispersion relation. The 
results for our three-channel fits are given by: 



,(6.10) 
2 

,(6.11) 



(0) 
(0) 



-0.195, 
-0.269, 



F (6.10) (()) 



,(6.11) 



(0) 



0.426 , 
0.424. 



(4.7) 



As can already be guessed from the similarity of the two- and three-channel form factors, 
F 3 (0) for the fits (6.10) and (6.11) comes out very close to the corresponding numbers 
in the two-channel case. The result for ^(0) should be compared with the tree-level 
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Figure 8: |F 3 (s)| for the fits (6.10K3) and (6.10K4) while varying Fi(A^). The dotted 
lines correspond to F 1 (A K7T ) = 1.21, the solid lines to 1.22, and the dashed lines to 1.23. 
The curves which are higher at the Kq(1A3Q) resonance always correspond to (6.10K4). 

expectation from U(3)lxU(3)_r xPT in the large- N c limit including resonances, F 2 XPT (0) = 
— /xr/(0)/3 = —0.109. Here for (6.11) the deviation amounts to more than a factor of 
two, and for (6.10) it is somewhat smaller. However, this difference is similar to that of 
-F 3 (0) before improving the high energy tail of the experimental data above 1.75 GeV, and 
hence, until such an improvement is performed for the three-channel case, one cannot draw 
a definite conclusion on ^(0). 

Like in the two-channel case, also for the three-channel fits (6.10) and (6.11) we inves- 
tigate the low-energy constraints Fi(Aktt), as well as the derivative of Fi(s) at zero -Di(O). 
The results are found as: 



Again, Fi(Ax n ) for three channels turns out very close to the two-channel case where 
also Fi(Axtt) has been found too large, so that the corresponding form factors cannot be 
considered physical. Nevertheless, the finding that for both fits Fi(Ax n ) is very close is 
reflected in the fact that the full form factors i*\(s) turn out very similar. The deriva- 
tive of Fi(s) at zero -Di(O) for (6.10) and (6.11), on the other hand, also is close to the 
corresponding numbers in the two-channel well as to the x?T expectation. 



if^A^) = 1.266, 
Fl 6A1 \A Kn ) = 1.264, 



Dj 6 - 10 )(0) = 0.927, 
£>{ 6 - n) (0) = 0.909. 



(4.8) 
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Figure 9: |i*i(s)| in the three-channel case. The thick solid and long-dashed lines corre- 
spond to the unitarised x?T fits of eqs. (6.10) and (6.11) of respectively. For comparison, 
as the thin lines we have also displayed the corresponding two-channel form factors for 
(6.10) and (6.11). 

5 Conclusions 

In this work, we have investigated the strangeness-changing scalar form factors for the 
Ktt, Kr] and Krj' systems up to 2 GeV. To this end, first theoretical expressions for 
the scalar form factors were derived in the framework of the chiral Lagrangian including 
resonances in the limit of a large number of colours. In a second step, analyticity and 
unitarity constraints were employed to write down a set of dispersion relations for the form 
factors. Taking into account previous results on S-wave Kir scattering fl[], the scalar Kir, 
Kr) and Krf form factors could be extracted from the solutions of the coupled-channel 
dispersion relations. 

The f] and rj' mesons can be described in the U(3)z,xU(3).r effective theory in the large- 
N c limit. In this framework, theoretical expressions for the tree-level scalar form factors 
have been derived in section 2. However, the couplings of the scalar meson resonances to 
the pseudoscalar mesons, q and c m , which appear in the resonance chiral Lagrangian are 
not very well know. Assuming the scalar form factors to vanish at infinity, constraints on 
the scalar couplings could be obtained which are compatible with present phenomenological 
knowledge on a and c m . In particular, if the chiral low-energy constants L5 and L$ are 
saturated by the scalar meson resonances, we deduced the relation L 5 = 2L 8 , independent 
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Figure 10: | -F2 ( -s ) | and |-Fs(s)| in the three-channel case. The solid and long-dashed lines 
correspond to the unitarised fits of eqs. (6.10) and (6.11) of Q respectively. The two 
lower lines correspond to |F 2 (s)|. For comparison, as the thin lines we have also displayed 
the corresponding two-channel form factors for (6.10) and (6.11). 

of the number of resonances. Further consequences of the theoretical constraints have 
already been investigated in our previous analysis of S-wave Kit scattering and in 
more detail in section 2.2. Furthermore, in the chiral framework it was found that the Krj 
form factor is strongly suppressed compared to the Kir and Krf form factors. 

Using unitarity and analyticity, we were able to write down a set of coupled disper- 
sion relations, which relate the scalar form factors to the S-wave, isospin 1/2 scattering 
amplitude. In the elastic, single-channel case, an analytical solution of the dispersion re- 
lation due to Omnes is known, which only depends on the elastic S-wave phase shift. In 
the coupled-channel case the set of integral equations had to be solved numerically. To 
perform this step, an efficient algorithm to perform the Cauchy principal value integrals 
was developed, which is described in appendix C.2. 

Using our results for different fits to the elastic Kn scattering data it was found 
that the scalar Ktt form factor Fi(s) is not well determined in the single channel case 
and that higher energy data, as well as inelastic channels have to be included for a better 
description. Since it was shown that the Kr\ channel is suppressed, we have concentrated 
our detailed investigation of the form factors on the two-channel case with Ktt and Kr( . 
In general, the two-channel dispersion relations admit two linearly independent solutions, 
and thus we need two integration constants to fix them unambiguously. 
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An obvious choice for the first constant is -Fi(O) which is very well known from next-to- 
leading order %PT. Different possible choices for the second constant could be ^(0), the 
derivative of Fi(s) at zero -Di(O), or ^(A^,,-) at the so-called Callan-Treiman point. We 
have used the last constraint because in our opinion it is known most precisely also from 
next-to-leading order %PT. With this information we then deduced F 3 (0) = 0.58 ± 0.04, 
compatible with the tree-level xPT expectation in the large- N c limit. Furthermore, -^3(0) 
is also in excellent agreement to the value from those fits with an improved high energy 
behaviour that gave rise to only one independent solution of the form factors. 

The three-channel case was used to corroborate that indeed the Krj form factor ^(s) 
turns out much smaller than Fi(s) and F 3 (s), and that the scalar Ktc and Krj' form factors 
were practically unchanged compared to the two-channel case. This gives us confidence 
that the more elaborate two-channel analysis is sufficient for our purpose of determining 
the strangeness-changing scalar spectral function since the Krj contribution is negligible. 

In a forthcoming work f52fl , our results for the scalar Kit and Krj' form factors will 
be utilised to explicitly calculate the corresponding strangeness-changing scalar spectral 
function up to 2 GeV. This spectral function is a key ingredient in the determination of the 
strange quark mass from QCD sum rules for the scalar correlator and it is also important 
in the corresponding analysis of the Cabibbo-suppressed hadronic r decay width. 
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Appendices 



A Loop functions 

For completeness, we tabulate here the one-loop function Jpq(s) [Q appearing in the next- 
to- leading order chiral form factors of eqs. (|2.7|) and fl2.8| ) : 



Jpq(s) 



1 



167T 2 

1 



dx In 



Mp — sx(l — x) — A P qx 
Ml - A PQ x 



(A.l) 



32vr 2 



A 



In 



' M. 



Q 



PQ, 



In 



[s + Xpq) 2 — A 



PQ 



Is -\ PQ ) -A 2 pq 



with 



S PQ = M 2 +M 2 ; A P q = M 2 -M 2 ; \ 2 PQ = 



s-{M P + M Q f s-{M P -M Q f 



B Three— channel dispersion relations 

Below we present the dispersion relations for the scalar form factors Fk(s) in the full 
three-channel case including Ktt, Ki] and Krj': 



F 2 (s) 



F*(s) 



+ 



+ 



7T 



Sthl 
oo 



71 



ds 



[s? -s-iO 



s th3 



— ds 

7T J 

Sthl 



(s' -s-iO 



s' -s- iO 



7T J (S ; — S — «0 

S«l3 



7T 



7T 



(/,S 



(s' - s - iO 

,a 3 (s')F 3 (s')tf 
is' -s- iO 



s th3 



s'Y , 1 J ds/ a 2 (s')F 2 (s')tl 2 (s f y 



ix J (s' — s — iO) 

Sth2 



7T 



s th2 



[s'-s- iO) 



7T 7 (s' — S — ZO) 

•5th 2 



(B.l) 
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C Numerical solution of the dispersion relation 



To illustrate our method for the numerical solution of the dispersion relation, below we shall 
discuss the single channel case. The coupled channel case can be treated in a completely 
analogous fashion. 



C.l Single channel case 

In the single channel case, the dispersion relation reads: 

ReF 1 (s) = vlds'-^i- = vjds'^^+vjds'/^, (C.l) 
J is' — s) J is' — s) J is' — s) 

s t h sth s cut 

/l(s ') = limbs') = -a^OFiOOCOO*, (C2) 

7T 7T 

and s t h = (Mk + M n ) 2 . For numerical purposes on the right-hand side we have split the 
integration range at an energy s cu t. Equation (|C.1|) can also be written as a once subtracted 



with 



dispersion relation: 

" jh(s') , sh(s') 



Sth s t h 



Assuming that fi{s) vanishes for s — * oo, such that the first integral converges, it can 
be identified with -Fi(O). In the numerical analysis, we have also investigated solving the 
subtracted dispersion relation imposing i<i(0) as a boundary condition. 

For the numerical integration, it is convenient to perform a change of variables such 
that the integration range becomes finite. We shall transform both energy intervals to the 
interval (0, 1). For the low-energy interval (sth, s cu t) a convenient transformation is: 



s' = s th + (scut - sth) x' 2 ; x' 



whereas for the high-energy region (s cu t, oo) we chose: 



v , (C4) 

\ {Scut — Sth) 



{l-bz') _ {s 1 - sth) (c , 

s = s th -r, ~pr , z — -r~i u V " V^-^) 

(1-2) {s'-bsth) 

The additional parameter b < 1 in the second transformation has been introduced to allow 
for increasing the convergence of the method or to check the stability of the solution. 
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Defining the additional variables x = \J (s — s t h)/(s cu t — Sth) and z = (s — s t h)/(s — bsth), 
the unsubtracted dispersion relation takes the form: 

whereas the subtracted dispersion relation is given by: 

Re Fl (x) = Fl (0) +Vfdx' + +v f dz , (l-bz)f^) 

with a = s cut / s t h- The additional factor of x' in the low-energy integral has the advantage 
of smoothing out the square-root singularity x = x' — 0. 

As a further ingredient, for the iteration of the dispersion relation, we require a rea- 
sonably accurate numerical integration routine for the Cauchy principal value integrals. 
Generally, the most accurate integration routines are based on the Gauss algorithm. In 
the next section, we shall thus develop a Gauss algorithm for the Cauchy kernel. 

C.2 Gauss quadrature for Cauchy kernel 

Gaussian quadrature algorithms can be designed such that the approximation 

N 



J f(y)w(y)dy « £ w sf(Vi) 

a J'=l 



is exact if f(x) is a polynomial. The task for finding such algorithms is finding the set of 
weights Wj and abscissas yj to accomplish this feat, given a weight function w(y). 

The starting point for developing the Gauss algorithm lies in finding a complete basis 
of functions which are orthogonal with respect to the weight function w(y) f55|. For the 



trivial weight w(y) = 1, on the interval (—1, 1), this set of functions can be chosen to be 
the Legendre polynomials P n {y) and the corresponding integration formula is the so-called 
Gauss-Legendre algorithm. 

Defining the scalar product of two functions / and g over the weight function w(y) as 

6 

(f\g) = J f(y)g(y)w(y)dy, (C.9) 

a 

we therefore seek functions u n (y) such that (u m \u n ) is zero if m ^ n. Choosing the Cauchy 
weight w(y) = l/(x — y), it can be shown that again on the interval (—1,1) a set of 
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orthogonal functions is given by 

u n (x,y) = P n (y) - 



Q, 



X 



Qn-l(x] 



P n -i{y), 



(CIO) 



where Q n (x) are the associated Legendre functions of the second kind [56]. The normali- 
sation is easily found to be: 



(u n \u n ) = 2Q n (x) 



Pjx) 



Qr 



[X 



Qn-l(x) 



Pn-l{x) 



2Q n (x) u(x,x) 



(C.ll) 



For a Gauss algorithm of order N, the set of abscissas is now given by the zeros yj of the 
function u^iy) in the interval (a,b), and the weights Wj can be calculated from a general 
formula, which can for example be found in the Numerical Recipes |55 |. In our particular 
case, however, we still have the complication that the weight function, and thus also the 
u n , also depend on the additional variable x. The problem simplifies considerably, if we 
only evaluate our integral at the zeros Xi of the function Q^(x). Then our set of orthogonal 
functions is again the Legendre polynomials, and from a straightforward calculation, the 
weights Wj(xi) are obtained to be 



Wj(Xi) 



2(1 - yj) 

N(xi - 



Q N -i(xj) P N (xj) 

P 2 N -i(yj) 



(C.12) 



Our final integration formula reads: 



f(y) 



dy 



N 



[Xi 



y) 



J2 w A x i)f(vj) 

3=1 



(C.13) 



It is a trivial matter to find the integration formula for a general interval (a, b) using 
linear transformations. The result for the integral at values of x different from X{ can 



be calculated with the help of standard interpolation algorithms ||55|| . A final comment 
concerns the practical calculation of the roots of Q^. Once the N roots of P/v o n the 



interval (—1,1) are computed with a standard algorithm [|55 |, the N + 1 roots of Qn are 
easily searched for, because the always interleave the roots of P^. 

Although the presented Gauss algorithm is rather simple, we were unable to find it in 
the literature. Thus, in appendix D we present our Fortran implementation of the algorithm 



which is based on the Gauss-Legendre algorithm given in the Numerical Recipes |55 . 
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D Fortran code for the Gauss routine 



C ! Small program to test the subroutine "cauleg" 

C ! Calculates the integral int( sin(y) ~2/(x-y) , y=0:3); 

C ! Last change: 12.1.2001 

C ! Matthias Jamin: m.jamin@thphys.uni-heidelberg.de 

Program CauLegTest 

Implicit Double Precision (A-z) 

Integer i, j , ngau 

Parameter (ngau =63, a = 0.D0, b = 3. DO) 
Dimension x(ngau+l) , y(ngau) , w(ngau+l ,ngau) 
! 

C ! Calculate zeros of the Legendre polynomials and Cauchy weights 

! 

Call cauleg (a, b,x,y,w, ngau) 
! 

C ! Calculate integral at the zeros x(i) of Q_n(x) 

! 

Open (unit=8, file=' cauleg. out ' ) 

Do 20 i=l,ngau+l 

! 

sum = 0.D0 
Do 10 j=l,ngau 
10 sum = sum + dsin(y ( j ) ) **2*w(i , j) 
sil = sum 



With the CERN library "mathlib" this is the exact result 



xx = x(i) 

si2 = ( dcos(2.D0*xx)*(dcosin(2.D0*(a-xx))-dcosin(2.D0*(b-xx))) 
+ dsin(2.D0*xx)*(dsinin(2.D0*(xx-a))-dsinin(2.D0*(xx-b))) 
+ dlog((b-xx)/(xx-a)) )/2.D0 

! 

20 Write (8,*) x(i) , sil, si2 
Close (8) 
! 

Return 
End 
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c 



C ! Gauss routine for integration with the Cauchy kernel 1/ (x-y) 

C ! based on the routine "gauleg" from the Numerical Recipes. 

C ! Uses legp(x,n) and legq(x,n) 

C ! Given the lower and upper limits of integration xl and x2, and 

C ! given n, the routine returns arrays x(l:n+l), y(l:n), w(l:n+l,l:n) 

C ! containing the points x where the integral should be evaluated, 

C ! as well as the abscissas and the weights of the Gauss-Legendre 

C ! n-point quadrature formula. x(l:n+l) and y(l:n) contain the roots 

C ! of the Legendre functions Q_n and P_n of order n, shifted to the 

C ! interval (xl,x2). Other values of x have to be calculated with 

C ! standard interpolation methods. 



Subroutine cauleg(xl,x2,x,y,w,n) 
Implicit Double Precision (A-z) 
Integer i , j , m, n 
Dimension x(n+l), y(n), w(n+l,n) 

Parameter (Pi = 3 . 141592653589793238D0, eps = l.D-15) 
! 

m = (n+l)/2 

xm = 0.5D0*(x2+xl) 

xl = 0.5D0*(x2-xl) 

! 

Do 20 i=l,m 

z = dcos(Pi*(i-0.25D0)/(n+0.5D0)) 
10 zl = z 

Call legp(pl,pp,z,n) 
z = zl-pl/pp 

If (dabs(z-zl) .gt.eps) Goto 10 
y(i) = xm-xl*z 
y(n-i+l) = xm+xl*z 
20 Continue 
! 

Do 40 i=l,m+l 

If (i.eq.l) Then 

z = -1.D0+1.D-9 
Else 

z = 0.5D0*(y(i)+y(i-l)-2.D0*xm)/xl 
End If 
30 zl = z 

Call legq(ql,qp,z,n) 
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z = zl-ql/qp 

If (dabs(z-zl) .gt.eps) Goto 30 
x(i) = xm+xl*z 
x(n-i+2) = xm-xl*z 
40 Continue 
! 

Do 50 i=l,m+l 

xi = (x(i)-xm)/xl 
Call legp(pl ,pp,xi ,n) 
Call legq(q2,qp,xi,n-l) 
Do 50 j=l,n 

yj = (y(j)-xm)/xl 

Call legp(p2,pp,yj ,n-l) 

w(i,j) = 2.D0/dble(n)*(l .D0-yj*yj) 

/(yj-xi)*pl*q2/p2/p2 
w(n-i+2,n-j+l) = -w(i,j) 
50 Continue 
! 

Return 
End 



! Calculates the Legendre polynomials P_n(x) and P'_n(x) 
! P_n(x) is in pi and P'_n(x) in pp. 

Subroutine legp(pl ,pp,x,n) 
Implicit Double Precision (A-z) 
Integer j , n 
! 

p2 = 1.D0 

pi = x 

Do 10 j=2,n 

p3 = p2 

p2 = pi 

pi = ((2.D0*j-l.D0)*x*p2-(j-l.D0)*p3)/j 
10 Continue 

pp = n*(x*pl-p2)/(x*x-l.D0) 
! 

Return 
End 



! Calculates the Legendre functions Q_n(x) and Q'_n(x). 
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C ! Q_n(x) is in ql and Q'_n(x) in qp. 

! 

Subroutine legq(ql ,qp,x,n) 
Implicit Double Precision (A-z) 
Integer j , n 
! 

q2 = 0.5D0*dlog(dabs((l.D0+x)/(l.D0-x))) 
ql = x*q2-l.D0 
Do 10 j=2,n 

q3 = q2 

q2 = ql 

ql = ((2.D0*j-l.D0)*x*q2-(j-l.D0)*q3)/j 
10 Continue 

qp = n*(x*ql-q2)/(x*x-l.D0) 
! 

Return 
End 

c 

I 
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